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ABSTRACT 


) 

A  multi-step  approach  to  rectify  radar  imap.es 
reometrical ly  and  radiometr ically  is  presented.  The 
-ethod  uses  digital  elevation  models  (DEMs)  and  an 
.image  simulator  to  establish  the  relation  between  radar 
-image  coordinates  and  ortho  image  locations.  The 
simulation  serves  to  create  a  synthetic  gray  value 
samage  similiar  in  its  geometry  to  the  real  radar  imagey 
—associating  with  each  image  pixel  a  backscatter  value 
dependent  on  the  incidence  angle  and  a  DEM  address. 
Correlation  between  simulated  and  real  image  results  in 
2  relation  between  radar  backscatter  values  and  DEM 
addresses  and  produces  an  ortho-image  in  an  additional 
step. 


Stereo-derived  and  map-derived  elevation  models 
were  used  to  demonstrate  the  rectification  process  for 
-air-  and  spaceborne  radar  images. 


1.0  INTRODUCTION 


Numerous  remote-sensing  applications  are  based  on 
precisely  relating  the  gray  values  of  an  aircraft  or 
satellite  image  to  the  corresponding  phenomenon  on  the 
ground,  to  gray  values  of  other  images,  or  to  symbolic 
descriptions  of  an  object  in  the  form  of  maps. 

One  concept  to  solve  this  correspondence  problem 
is  by  physically  generating  a  derived  image  that  is 
geometrically  corrected  and  therefore  has  the  geometric 
properties  of  a  map.  The  result  can  be  denoted  by 
"ortho-images1'  or  by  "geo-coded  imagery".  The  first 
expression  results  from  standard  photogrammetrlc 
terminology  where  sinoe  about  1955  the  ortho-photo  has 
become  a  routine  product  (Manual  of  Photogrammetry , 
19411).  "Geo-coding"  is  a  more  recent  terminology  that 
was  Independently  developed  in  the  remote-sensing  and 
image  processing  fields  (Guertin,  Shaw,  1981). 

The  generation  of  digital  geo-coded  ortho-images 
requires  that  given  raw  input  images  are  resampled  In  a 
digital  image  processing  system.  Resampling  must 
result  in  a  new  output  image  so  that  the  image  gray 
value  of  each  imaged  object  is  presented  in  its  proper 
position.  Essentially  therefore  one  has  to  shift  input 
image  positions  x,  y  into  output  positions  X,  Y  such 
that 

-  effects  of  the  sensor  position  and  attitude  during 
image  acquisition; 

-  effects  of  perturbations  of  the  internal  sensor 
operations ; 

-  effects  of  the  imaging  characteristics  or  type  of 
geometric  projection; 

-  effects  of  the  object  shape 

are  all  eliminated  from  the  image  geometry.  The 
resulting  output  image  has  a  specified  geometry  that  is 
independent  of  the  specifics  of  sensor,  and  image 
acquisition. 
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This  report  describes  one  procedure  to  take  a 
digital  radar  image,  recordings  of  the  position  and 
attitude  of  the  sensor  platform,  and  a  description  of 
terrain  shape  in  the  form  of  a  digital  elevation  model 
(DEM),  to  obtain  a  digital  ortho-radar  image.  The 
approach  is  based  on  a  pixelwlse  re-arrangement  of  the 
input  data;  therefore  one  can  refer  to  the  procedure 
as  "differential  rectification". 

As  a  by-product  of  geometric  rectification  one  can 
also  manipulate  the  image  gray  values  themselves  in  a 
step  of  "radiometric  rectification",  thereby 
eliminating  image  gray  value  effects  of  topographic 
terrain  slopes. 

The  report  will  first  review  the  current 
understanding  of  radar  rectification  and  single  image 
radargrammetry .  It  will  then  discuss  the  individual 
steps  of  the  rectification  procedure  developed  in  the 
current  context,  and  finally  demonstrate  results  using 
aircraft  and  satellite  images,  both  with  map-derived 
and  with  stereo-radar  terrain  elevation  data. 


REVIEW  OF  SINGLE  IMAGE  RADARGRAMMETRY 
RECTIFICATION 


2.1  Inage  Coordinates 


Radargrannetry  requires  first  of  all  that  the 
radar  measurements,  tine  t  and  slant  range  r,  are 
properly  extracted  from  the  radar  recording.  This 
inplies  some  transformation  of  an  image  coordinate 
pair,  x,  y  or  pixel  line  and  sample  numbers,  i»j»  into 
r  and  t.  The  basic  concept  is  illustrated  in  Figure 
2.1. 

We  measure  for  each  image  point  P  its  x^y^ 
comparator  coordinates  or,  if  we  deal  with  a  digital 
image,  its  corresponding  line  and  sample  numbers. 
These  need  to  be  converted  to  the  time,  t,  when  object 
P  was  imaged  (sensor  events  time),  and  range,  r,  to  P. 
Numerous  methods  exist  to  relate  xj(,y|(  to  r,t. 
Generally  we  may  therefore  Just  define  an  operator  Ft 


r*=F,(x*.y*) 
t  =  F,(xk,y„) 


(2.1) 


X  \>: 


Neaavred 
Image  Petal 


FIGURE  2.1i  Measuring  coordinates  in  a  side- 

looking  radar  image;  x^y^  are  compa< 
rator  coordinates,  r,t  are  range  and 
tine  coordinates. 
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In  certain  types  of  radar  images,  particularly 
those  designed  for  good  geometry,  one  may  find 
auxiliary  reference  marks,  such  as  those  of  Figure  2.1, 
denoted  by  Rj,  R2  etc.  Depending  on  the  type  of 
marks  we  can  define  the  form  of  Equ.  (2.1)  and  its 
parameters.  As  an  example,  a  single  row  of  so-called 
time-marks  may  enable  one  to  define  a  straight  line  to 
be  denoted  as  "image  time  coordinate"  axis,  x. 

Numerous  methods  could  be  defined  to  actually 
implement  operators  Fr,  F  t .  Factors  to  be 
considered  are  the  type  of  geometric  referencing,  the 
type  of  radar  record  in  a  slant  or  ground  range 
presentation,  the  individual  digital  or  electronic 
pre-processing  of  the  image  data  etc. 

In  most  cases,  radargrammetry  begins  with  the 
SLR-image,  not  with  the  raw  electronic  signals  acquired 
by  the  antenna. 

However,  a  trend  in  digital  radar  systems  is 
developing  to  be  beyond  the  radar  record  and  to  look 
into  raw,  unprocessed  signals.  First  indications  are 
visible  in  studies  by  Kratky  (1979)*  Brown  et  al. 
(1981)  and  Curlander  et  al.  (1981).  These  new 
approaches  rely  on  concepts  such  as  "sensor  events 
time",  l.e.  the  exact  time  when  a  raw  pulse  was 
transmitted  from  and  received  by  the  radar  antenna,  and 
how  the  raw  radar  echo  is  converted  into  image  points. 

Data  from  the  SEASAT-SAR  experiment  can  serve  to 
illustrate  the  problem  of  extracting  t,  r  from  the 
image. 

For  a  long  time  the  digitally  correlated 
SEASAT-SAR  images  did  not  show  any  time  or  range 
reference  marks.  The  user  was  not  provided  with  means 
to  enable  him  to  relate,  x,  y  to  t,  r  coordinates. 
However,  Curlander  (1981)  recently  established  a 
relationship  that  exists  between,  on  one  hand,  the 
sensor  events  time  t  and  slant  range  r  and  on  the  other 
hand,  the  image  pixel  line  (azimuth)  number  1  and 
sample  (range)  number  J.  This  relationship  is  for 
digitally  correlated  SEASAT-SAR  rather  complex.  This 
is  illustrated  by  the  formulas  that  were  presented  by 
Curlander  (1981),  who  vent  back  to  the  raw  signal  data 
to  define  r,  t  for  each  pixel. 
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2 . 2  Projection  Equation 

Numerous  authors  have  formulated  rigorous  radar 
projection  equations,  similar  to  the  so-called  rigorous 
photogrammetrlc  equations  for  photography.  The  intent 
of  these  formulations  is  to  model  the  typical  range 
projection.  Object  points  are  projected  along 
concentric  circles  onto  a  projection  plane  (Figure 
2.2).  Authors  who  have  contributed  here  are  Rinner 
(1948),  Konecny  and  Derenyi  (1966),  Rosenfield  (1968), 
Akowetzki  (  1968),  Derenyi  (  1970),  Grade  et  al. 
(1970),  Leberl  (1970),  Hockeborn  (1971),  Norvelle 
(1971),  Geier  (1972),  Greve  and  Cooney  (1974), 
DBA-Systems  (1974)  and  many  authors  since.  A  review 
was  compiled  by  Leberl  (1976b). 
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FIGURE  2.2:  Range  projection.  Each  image  point 

defines  a  circle  as  the  locus  of  the 
corresponding  object  point. 


For  such  a  formulation  we  need  to  define  the 
sensor  position,  s(t)  e  (sx(t),  Sy(t),  s2(t))T, 
and  antenna  attitude  with  $(tj,  u(t),  <(t).  In 

synthetic  aperture  radar,  tue  synthetic  antenna’s 
attitude  is  merely  the  first  derivative  of  position, 
namely  the  velocity  vector  s(t)  =  (sx(t),  sy(t), 
Az(t)).  Both  s ( t )  and  s(t)  are  functions  of  time  t. 
The  formulation  also  requires  a  definition  of  objeet- 
and  sensor  coordinate  systems  according  to  Figure  2.3 
with  unit  vectors  e1f  e2,  e3  and  u,  v,  k, 
respectively. 

He  then  find  for  an  object  point  P  its  position 
vector  £: 

P=P,£.+P,£:  +  P«£s  ro 


and  in  the  sensor  system 


5 


4 


and 


FIGURE  2.3:  Definition  of  object  space 
sensor  coordinate  systems. 


Therefore, 


p=s+r 


\ 

(2.4) 


It  is  straight  forward  to  see  from  Figure  2.4  that 


u,=r  sin  t 

t>J,=r(sin2  fl-sin2  t),,j 
wr  =  -  r  cos  A 


(2.5) 


I 


Auxiliary  angle  *  is  a  system  constant  called 
squint-angle  by  most  authors  (Rosenfield,  1968; 
Hockeborn,  1971)  and  defines  the  conical  shape  of  the 
radar  pulse. 
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Introducing  a  new  vector  £  =  (up,  vp,  wp)  in 
the  sensor  coordinate  system  and  a  rotation  matrix  A  to 
rotate  the  sensor  system  u ,  v ,  w  into  the  object  system 
»  ®2»  £3»  w®  find  in  matrix  notation: 


FIGURE  2. A:  Relating  squint  angle  t,  slant  range 
r  and  elevation  angle  U  to  sensor 
coordinates  u,  v,  w. 

P=s+A-p  (2.6) 

where 

MO  «,(»)  MO" 

p,(t)  vfl)  p,(f)  (2.7) 

*M0  w,(0  w^f) 

Rotation  matrix  A  can  contain  the  classical  orientation 
angles  ♦  ,  «  and  "k  if  we  deal  with  a  real  antenna.  With 
synthetic  antennas,  we  obtain: 
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u  =  (u„,  ur  u,)  =  s/|s| 
r  =  (i’x,  vr  v ,)  =  (sx  j)/|s  x  sf 
w  =  (wx,  wr  w,)=(ux  v)/l u  X  tf| 


(2.8) 


The  unit  vector  u  extends  along  the  synthetic  antenna, 
thus  in  the  direction  of  the  velocity  vector  s.  The 
unit  vector  v  should  be  normal  to  u  and  to  the  nadir 
direction  of  the  antenna.  For  satellite  sensor,  the 
nadir  direction  is  easily  seen  to  be  the  direction  of 
the  sensor  position  vector  in  a  planetocentric  system. 
With  aircraft  radar,  one  has  to  specify  here  explicitly 
that  a  geocentric  system  is  assumed  also  for  the  sensor 
position  vector  s  and  the  object  coordinate  system 
defined  by  et,  e2,  e3. 

Equations  (2.5)  and  (2.6)  relate  object  and  image, 
where  r  and  time  t  derive  from  the  image  coordinates. 
We  now  see  that  -  contrary  to  classical  photogrammetry 
-  we  do  not  algebraically  relate  image  and  object 
coordinates,  but  the  image  only  defines  time,  which  in 
turn  defines  the  parameters  s  and  A  of  Eq.  (2.6). 

If  the  imaging  arrangement,  and  thus  s(t),  A(t), 

are  known  (measurement  of  navigation  data),  then  we 
have  for  each  point  3  equations  (Eq.  2.6)  with  4 
unknowns  i.e.  fl  ,  px,  Py  and  ps.  It  may  be 
preferable  to  work  with  2  equations  and  3  unknowns. 
This  is  achieved  by  elimination  of  the  unknown 
elevation  angle  ft  (angle  between  nadir  and  line  of 
sight),  thus  giving  the  following  two  equations: 


Ip  -  s| = r 


u  (p-s)  =  sin  T|u||p-s|  =  sin  tt 


(2.9) 


(2.10) 


where  we  now  have  a  simple  geometric  description  of  the 
radar  projection  case.  The  geometric  locus  of  a  point 
is  the  intersection  of  range-sphere  (Eq.  2.9)  with  a 
cone  (Eq.  2.10).  This  is  the  circular  projection  line 
mentioned  before. 
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From  Figure  2.9  one  finds: 

(s,-p)  5,  =  |s(-£|  Isjeosa, 


where  is  the  angle  between  the  linear  orbit 

element  i  and  the  vector  from  a±  to  .  It  relates 
to  entity  egi 

cosa  =  e,/|s,-p| 


where  eg  results  from: 

(£,-£)•£,  =  Is,!  *,,  e, =s4  (s1-p)/|si| 


With  this  we  get  orbit  position  sp: 
s,=sl+i1-e#/|il| 


(2.25) 


Position  Sp  then  leads  to  slant  range  rp  and  time 
tp  of  imaging: 

r,=(p-J,)  (2.26) 


l  =  f  ( + Is,  -  s,|  •  (» i + ,  - 1,  )/|sJ  (2.27) 

Actual  image  coordinates  xp,  yp  are  obtained 
through  an  inversion  of  the  equations  of  inner 
orientation : 

x,  = x,  +  (t  - (x, + ,  -  x,)/{tt  +  1  -  f,) 

y, =r/f-c 


Here,  x*f  Xi+j  are  image  x-coordlnates  of  time 
references  t*, 

These  computed  image  coordinates  can  be  compared 
with  actual  measurements  of  x,  y.  Discrepancies  are 
the  consequences  of  errors  of  inner  and  exterior 
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require  correction  with  the  help  of  ground  control 
points. 

Figure  2.9  illustrates  the  solution  to  the 
computation  of  the  image  coordinates  x,  y  for  a  given 
object  point  P(px»  Py,  pz)  as  it  was  published  by 
Leberl»  Fuchs  and  Ford  (1981).  The  exterior 
orientation  is  known  in  the  form  of  position  vectors 
*i »  i=1»...,m  at  times  t^,  1=1, Also  given 

is  thus  the  velocity  vector  s^  s  s^^.^  -  s <  ,  or 

the  equivalent  attitude  vector  describing  4>(t),  <(t). 


times  t i f  t2,...|  at  time  tn,  the  orientation 
parameter  is  considered  to  be  a  linear  extrapolation 
from  m  previous  times  plus  a  random  addition,  e.  Ebner 
and  Hoessler  could  show  that  a  numerically  stable  and 
comparatively  manageable  solution  can  be  formulated.  A 
detailed  description  is  beyond  the  scope  of  this  paper. 
An  application  to  radar  must  still  be  investigated. 


2.U.3  Hybrid  Methods  - 

A  combination  of  parametric  and  interpolative 
methods  is  feasible  and  common  in  photogrammetry .  A 
mathematically  rigorous  model  of  the  imaging  process* 
such  as  used  in  the  parametric  case*  may  not  be  able  to 
describe  the  entire  chain  of  physical  events  leading  to 
the  image.  Unknown  elements  or  effects  must  be 
considered.  This  can  be  done  by  Introducing  Ax, 
Ay-corrections  to  the  image  coordinates}  these 
corrections  can  be  represented  by  polynomials  or  other 
interpolative  methods. 

In  photogrammetry  this  approach  is  denoted  by 
"self-calibration" y  the  concept  is  related  to  the 
method  of  so-called  "collocation"  (Moritz,  1973).  An 
elementary  implementation  is  by  solving  first  a 

parametric  formulation.  This  will  be  done  in  a 

least-squares  method  that  leaves  residuals  at  the 
observed  image  points.  The  residuals  are  then  subject 
to  an  interpolative  computation.  This  is  meaningful 
only  when  the  residuals  form  a  systematic  pattern.  A 

new  point  P  (xp,  yp)  is  thus  converted  to  object 
coordinates  by  first  adding  Axpy  Ayp  to  xpy  ypy 
and  by  then  applying  the  radar  equations. 


2.5  Projection  Of  A  Ground  Control  Point  Into  A  Radar 
Image  *" 

In  the  event  that  the  exterior  orientation  is 
measured,  e.g.  with  the  help  of  inertial  navigation, 
one  must  expect  errors  of  measurement  to  exist.  These 
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Instead  of  solving  d4> ,  d* ,  dsx,  ds  and  ds2  one 
solves  for  unknown  polynomials  coefficients  sq, 

Eq.  (2.21)  is  thus  rewritten: 

dx~  -{c0+c,t+c2fJ  +  •••} 

+P/{b0+b,t+b2t2  +  -  ••} 

-H{a0+a,f  +  a2fJ-f 

dy=-^ {d0+dtt+d2t2  +  •••}  i  (2'2 


+--{«o+es*+<jfJ  +  ---} 


Not  all  coefficients  can  be  determined  separately. 
The  analysis  of  this  approach  is  beyond  the  scope  of 
this  paper.  It  may,  however,  be  relevant  to  point  out 
that  instead  of  single  high  order  polynomials  it  is 
preferable  to  use  piecewise  polynomials  of  low  order. 
Baker  and  Mikhail  (1974)  have  studied  these  approaches, 
however,  for  aircraft  scanner  imagery,  not  for  radar. 

Dowideit  (1977)  studied  the  method  of  modelling 
the  exterior  orientation  parameters  by  Fourier  series, 
correctly  arguing  that  the  elements  of  the  exterior 
orientation  have  a  periodical  behavior  with  several 
dominating  frequencies  (Leberl,  1972a).  Again,  one  has 
to  solve  for  unknown  coefficients  of  the  Fourier 
series.  The  main  problem  is  instability  of  the 
resulting  equation  system  due  to  correlation  among 
coefficients . 

A  somewhat  different  approach  has  been  proposed 
for  dynamic  scanner  images  by  Ebner  and  Hoessler 
(1978).  The  elements  of  exterior  orientation  are 
considered  to  form  time-series  appropriately  modelled 
by  Gausa-Markof f-chains  of  order  m: 

-  ! 

0O.+i)=  I  +  i  j  (2.24) 

t*i  I 

i 

Thus  the  exterior  orientation  is  described  at  discrete 
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<ty*=sinn(- do)wp- dsf)-cos  (Hdu)-vf-ds,) 


iyms  +do)  sin  fi-rcos  O-sin  tlds, 
-  cos  0  du)  -r  sin  ft  +  cos  ft  <fs, 
</>>=  -sin  O  ds,+ cos  (I  ds, 
dx  =  -  r  sin  0  die  -  r  •  cos  Cldd>-  dsM 


(2.20) 


Vlth  flying  height  H  r  p  .  cos  0  ,  or  with  object 
coordinates  py  *  r  .  sin  0  *  this  can  be  rewritten: 


<fx*=  -dsM+p,-dK-H-d$ 
dy =  —pjr  dSj+H/rds, 


He  see  again  that  we  have  only  five  elements  of 
exterior  orientation.  The  sixth  element*  dw,  has  no 
effect  on  image  geometry. 

For  each  lmage-x~coordinate*  a  complete  set  of 
five  orientation  elements  must  be  solved.  This  is  not 
possible  with  a  co-linear  arrangement  of  control 
points*  since  clearly  the  effects  of  dsx  and  d<t> 
cannot  be  separated.  We  can  solve  either  (dsx* 
dsy,  dsx,  die)  or  <d$,  dsy*  dsx,  die).  The  fifth 
element  must  be  known.  The  required  number  of  control 
points  to  compute  the  four  elements  is  two. 

It  is  unrealistic*  however*  to  assume  that  for 
each  image  x-line  one  can  have  two  control  points. 
Also  this  is  not  necessary  because  of  correlation  that 
exists  between  the  exterior  orientation  values  at  two 
different  times  tj  and  tg» 

Authors  have  in  the  past  proposed  various 
mathematical  models  for  the  exterior  orientation  as  a 
function  of  time*  e.g.  DBA-Systems  (1974)  employed 
polynomials  such  that: 

d^(i)-fl0+ai*+fljfJ+"*  ! 


<k»(0"/o+/i*+/*r,+ '  *  * 
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(2.22) 


Differentiation  of  Eq.  (2.16)  leads  tot 


dp-dA  * (p  —  j)— A  ''Js  ^ 2  17) 

where  ds  =  (dsXl  dsy,  dst).  Since  we  assume 
♦  siiisit  0,  we  get  A'1  to  be  the  unit  matrix,  E. 

The  differential  form  d A” 1  .  (j>  -  s)  is  well 

known  in  photogrammetry  (Rinner  and  Burkhardt,  1972): 


dA~*  •(£-*)= 


0 

dx 


—dx 


0 


d<f> 

—do) 


A~‘(p-s) 


—d<t>  do)  0 


(2.18) 


Therefore  we  get  with  A“'  .  (£  -  s)  s  (up,  vp, 

wp)r, 


du. 

—  dKVp+d<fr-W0-dSM~ 

dv. 

- 

dK-Up-do)-Wp-dsr 

d»,_ 

-d4>-up+do)-vp-dst 

Errors  of  image  x , y-coordinates  result  from  dup, 
dvp  and  dwp  as  follows: 

dx=*duf 

But  Therefore: 


cx»s(ldwp 

Let  us  now  simplify  further  with  squint  angle  TsO. 
Then  ups0,  vp*r  .  sinD  and  wp«-r  .  coafl  .  Thus 


this  polynomial  need  to  be  computed: 

A  X,=a0+a,Xy+  •  •  • 

A  Y,=b0  +  blY,'+  - 


since  Xp,  Yp  s  0. 

This  computation  and  control  point  selection  is 
repeated  for  each  new  point.  In  the  case  of  a  first 
order  polynomial  one  also  has  the  expression  "moving 
plane”  interpolation.  A  zero-order  polynomial  leads  to 
a  simple  weighted  average  interpolation. 


2.11.2  Parametric  Models  - 

The  exterior  orientation  of  the  sensor  is  a 
function  of  time  and  therefore  continuously  changing. 
As  a  result  each  image-x-coordinate  relates  to  a 
different  exterior  orientation.  It  is  obvious  that 
there  cannot  be  rigorous  reconstruction  of  the  exterior 
orientation  from  control  points  for  all  times  t:  the 
number  of  unknown  exterior  orientation  elements  is 
infinite.  Even  if  we  had  several  control  points  imaged 
at  one  given  time  t  and  thus  along  one  line  in  the 
image-y-direction,  complete  computation  of  the  exterior 
orientation  is  not  possible. 

To  illustrate  the  problem  let  us  assume  that  the 
exterior  orientation  is  nearly  perfect  and  that  the 
flight  only  slightly  deviates  from  a  straight  line 
along  the  object  e^-coordinate  axis.  Also  t,  u,  k, 
the  elements  of  the  sensor  attitude,  are  nearly  zero. 

These  small  deviations  d$,  du,  die,  dsx,  dsy 
and  dsz  of  the  sensor  from  a  predictable  exterior 
orientation  can  be  related  to  image  coordinate  errors, 
dr,  dt,  or  dx,  dy.  We  start  from  Eqs.  (2.5)  and 
(2.6): 

i 

P-A-'ip-s)  |  (2.16) 
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FIGURE  2.7*  Moving  average  interpolation  with 

square  windows  to  collect  n  reference 
points. 


FIGURE  2.8:  Selection  of  control  points  in 
quadrants,  in  this  example  with 
2  points  per  quadrant. 

With  these  selected  control  point  sets  a 
polynomial  of  low  order  is  now  computed  (1st  and  2nd 
order)  whereby  each  control  point  obtains  a  weight  as  a 
function  of  its  distance  from  the  new  point.  If  all 
coordinates  have  been  reduced  to  an  origin  in  the  new 
point,  then  only  the  first  coefficients  an,  bn  of 


The  coefficients  ao,...,an,  bQ,...,bn  are 

computed  from  the  known  control  point  values  AX*, 
AY* . 

Various  authors  have  evaluated  different  types  of 
polynomials  for  radar  image  correction  (Leberlt  1970; 
Derenyi,  1974a;  Hirsch  and  van  Kullenburg,  1976; 
Shakine  and  LeToan,  1978).  A  common  conclusion  is 
difficult  to  draw,  since  different  data  from  different 
sensors  were  employed.  However,  one  may  summarize  that 
the  use  of  polynomials  is  generally  not  advisable  due 
to  their  predictably  unfavourable  behavior  in  the  case 
of  Irregular  control  point  distribution,  and  in  the 
case  of  extrapolation.  Polynomials  are  only 
recommended  if  they  are  of  low  order  (1st  or  second) 
and  if  there  is  an  even  control  point  distribution. 


(b)  Moving  average 

The  polynomials  of  Eq.  (2.15)  are  a  global, 
closed  representation  of  image  deformation.  There  are 
also  polnt~wise  methods  of  describing  this  deformation. 
We  start  from  the  known  set  or  AXj,  AYi-values  in 
control  points,  and  need  to  compute  the  AXp,  AYp  in 
a  new  point  Xp,  Yp.  The  moving  average  technique 
assembles  for  each  point  P  a  set  of  known  control 
points  according  to  some  criterion.  This  may  for 
example  be  the  n  closest  control  points  (Figure  2.7), 
or  it  may  be  m  control  points  in  each  of  four  quadrants 
(Figure  2.8). 


(a)  Polynomials 


He  may  now,  as  one  example,  attempt  a  scale  fit: 


(2.13) 


FIGURE  2.6:  Interpolative  method  of  relating  image 
and  map  points;  a  preliminary  trans¬ 
formation  leads  to  residuals  AX,  AY. 
These  define  a  correction  surface  that 
can  be  described  by  isolines,  separa¬ 
tely  for  AX  and  for  AY. 


Since  this  is  solved  with  more  control  points  than 
needed  one  obtains  residuals  AX,  AY  or,  if  Eq.  (2.13) 
is  applied  Inversely,  Ax,  Ay.  For  a  new  point, 
xp»Yp  we  find  thus  first  Xp,  Yp  and  we  obtain 
the  final  map  coordinates  by  adding  AXp,  AYp: 


+ A*, 

r,=Y;+AY, 


(2.14) 


where  AXp,  AYp  result  from: 

AX,«a0+a.*,+a2i;+a,X,i;+fl4  Y,J  +  •  •  • 

AY'-bo+^XsbiY'+btXJ'+btX,2*-  •  •  |  (2.15) 


(a)  Polynomials 


He  may  now, 


as  one  example,  attempt  a  scale 


fit: 

(2.13) 


FIGURE  2.6:  Interpolative  method  of  relating  image 
and  map  points;  a  preliminary  trans¬ 
formation  leads  to  residuals  AX,  AT. 
These  define  a  correction  surface  that 
oan  be  described  by  isolines,  separa¬ 
tely  for  AX  and  for  AY. 


Since  this  is  solved  with  more  control  points  than 
needed  one  obtains  residuals  AX,  AY  or,  if  Eq.  (2.13) 
is  applied  inversely,  Ax,  Ay.  For  a  new  point, 
xp»Yp  we  find  thus  first  Xp,  Yp  and  we  obtain 
the  final  map  coordinates  by  adding  AXp,  AYp: 


y,=  y;+Ay, 


(  2 . 1  A  ) 


where  AXp,  AYp  result  from: 


(2.15) 


A'-y-V- 


I 

Ay,-60+b1if,+6Jy,+6,ir,y,+64K/+ •••  | 


With  a  given  radar  image  and  collateral  imaging 
data  we  have  several  options  to  relate  the  image  to  the 
map  or  object.  We  generally  will  have  one  unknown  or 
at  least  only  erroneous  exterior  orientation.  This 
leads  to  the  necessity  of  using  control  points  or 
more  generally  -  control  information  to  relate  image 
and  object  space.  This  has  been  attempted  by  many 
authors.  The  efforts  can  be  classified  into: 

a)  lnterpolative  models  employing  so-called  "rubber- 
sheeting"  or  warping  algorithms  to  make  the  image 
fit  to  ground  control; 

b)  parametric  models  employing  projection  equations 
and  computing  unknown  parameters  from  control 
point  information;  and 

o)  hybrid  combination  of  a)  and  b)*  where  the  pro¬ 
jection  equations  are  used  to  relate  image  and 
objecty  but  where  an  lnterpolative  technique  is 
applied  to  represent  unknown  deviations  of 
reality  from  the  assumed  model. 


2.4.1  lnterpolative  Methods  - 

A  first  approach  to  the  study  of  image  geometry  is 
to  consider  the  image  itself  as  a  map  substitute  and  to 
analyse  its  differences  from  a  map  of  known  geometric 
properties.  Differences  Ax^,  AyA  between  image 
xi»  Yi  and  map  Xi*  Yi  coordinates  can  be 
measured  at  control  points  i.  Some  interpolation 
function  is  then  used  to  predict  Axp*  Ay- 
deformations  of  the  image  at  arbitrary  image  locations 
p.  The  same  methods  can  be  employed  as  those  used  in 
other  photogrammetric  interpolation  tasks  (Leberly 
1975).  To  give  but  two  examples*  let  us  consider 
Figure  2.6.  The  points  Pi*  i«1y...*n  are  known  with 
their  image  coordinates  x^y  yi*  and  map  coordinates 
Xif  Yi. 


unknown  coefficients  and  the  question  is  how  to  obtain 
estimates  of  these  coefficients.  As  will  be  shown 
later,  these  can  derive  from  known  object  point 
coordinates.  Error  mod-la  generally  imply  a  polynomial 
form  of  the  image  coordinate  deformation  Ar,  At. 
Grade  et  al.  (  1970),  Leberl  (1971),  DBA-Systems 
(197^)  and  others  have  applied  the  following 
polynomials : 

A/^ajr+flji^+ajr3  +  •  •  • 

A/=A0  +  Vf+A2.f*  +  A,.f  (2.12) 


Little  research  exists  on  the  various  alternatives 
available  to  describe  Ar,  At-errors  in  an  optimum  form. 
At  this  point  one  can  therefore  only  direct  attention 
to  the  fact  that  these  errors  must  be  somehow 
considered  in  radargrammetry  and  image  rectification. 


* 


FIGURE  2.5:  Radar  interior  orientation  and 
image  deformations  Ar,  At. 


In  Eq.  (2.10)  the  sensor  coordinate  axes  are  all 
described  by  the  geometric  position  and  velocity  vector 
of  the  antenna.  This  applies  in  particular  for  SAR, 
since  the  antenna  attitude  is  defined  by  the  velocity 
vector,  (the  "synthetic"  antenna  extends  in  the 

direction  of  the  physical  antenna's  motion).  For  real 
aperture  radar,  the  real  antenna's  attitude  itself  is 
to  be  used  to  define  u,  v  and  k.  Therefore  angles  4 
and  <  must  be  entered  in  Eq.  (2.7);  u  is 
insignificant  and  can  be  set  at  zero.  This  is  obvious 
from  Eqs.  (2.9)  and  (2.10)  which  do  not  contain  an 
angle  u  since  the  components  of  vector  u  in  the 
rotation  matrix  A,  namely  ux,  uy  and  uz  do  not 
contain  10  .  The  same  fact  is  also  obvious  from 
geometry:  the  projection  circle  is  not  affected  by  a 

sensor  roll,  u>. 

In  aircraft  radar,  t  normally  is  zero,  so  the  cone 
of  Eq.  (2.10)  degenerates  to  a  scanning  plane: 

u(p-s)= 0  (2.11) 

The  range  sphere  has  an  obvious  physical  explanation: 
the  echo-time  measurement  of  radar  clearly  associates 
with  each  image  point  a  slant  range  r,  and  therefore  a 
sphere,  as  one  geometric  locus  for  the  corresponding 
object. 


2.3  Inner  Orientation 

The  parameters  of  the  relationship  between  image 
coordinates  x,  y  and  the  object  space  entities  r,  t  are 
denoted  as  parameters  of  interior  orientation. 
Equation  (2.1)  was  a  mathematical  representation  of 
this  interior  orientation.  The  concept  derives  from 
photogrammetric  cameras,  where  each  image  point  defines 
a  projection  ray  attached  to  the  camera.  The  exterior 
orientation  relates  the  projection  ray  to  the  object 
point. 

However,  Eq.  (2.1)  represents  merely  a  model  of 
reality.  Errors  of  this  model  must  be  anticipated  as 
shown  in  Figure  2.5.  Of  course,  any  error  model  has 


orientation  x 

Ax=x-x 
Ay ’ey— y 

The  Ax,  Ay. discrepancies  in  control  points  can  be  used 
to  compute  correction  functions  and  then  to  compensate 
for  these  errors  in  new  points.  The  correction 
functions  can  be  of  the  types  as  described  in  Section 
2.4.1. 


2.6  Rectification 

In  orthophoto-generatlon  by  both  digital  image 
processing  and  photogrammetric  computer  assisted 
orthoprojection,  one  needs  to  define  a  sufficiently 
dense,  geometrically  correct  grid  of  anchor  points 
according  to  Figure  2.10.  These  anchor  points  are 
generated  using  rigorous  radargrammetry .  All  image 
details  inside  a  quadrangular  grid  mesh  are  transformed 
with  a  simpler  and  thus  faster  linear  method. 

y _ _  y 

,  i  i  i 

tzrf 

1  > 

i —  r 

X 

FIGURE  2.10:  Anchor  points  in  the  deformed  image 
(left)  and  orthophoto  (right)  (from 
Kraus,  1976). 

The  process  is  based  on  a  raster  DHM  at  orthophoto 
anchor  points.  These  terrain  points  are  resected  onto 
the  image  according  to  Eqs.  (2.25)  to  (2.27),  so  that 
a  situation  as  shown  in  Figure  2.10  results.  Based  on 
these  image  coordinates  one  of  the  standard  digital  or 
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analog  photogrammetr ic  rectification  routines  can  be 
applied . 

Rectified  aircraft  radar  images  were  obtained  from 
a  study  by  Leberl  and  Fuchs  (1978),  and  by  Leberl, 
Fuchs  and  Ford  (1981)  with  a  detailed  error  analysis. 
Its  discussion  here  would  go  beyond  the  scope  of  this 
report. 

Digital  rectification  has  been  done  so  far  in  the 
context  of  the  SEASAT-mission .  Curlander  et  al. 
(1981)  have  used  digital  SAR-images  from  SEASAT 
together  with  orbital  data,  and  Naraghl  et  al.  (1981) 
matched  the  SAR-data  with  a  digital  terrain  model. 

Instead  of  actually  creating  a  geometrically 
corrected  radar  image  one  can  think  of  converting  only 
certain  interpreted  image  details  to  a  map  projection. 
If  a  DHH  exists  with  the  image  and  the  image  is  related 
to  the  ground  height  through  a  computation  of  exterior 
orientation,  then  one  can  trace  the  image  lines  of 
interest  with  a  digitizer;  the  image  x,  y-coordinates 
are  transformed  to  X,  Y  knowing  Z.  This  concept  was 
realized  by  Greve  and  Cooney  (1974)  with  radar,  and  is 
a  method  of  photogrammetry  which  has  also  been  proposed 
by  other  authors  (e.g.  Makarovic,  1973). 

The  idea  has  not  yet  had  responses  of  any 
significance,  either  with  metric  photography  or  with 
radar  or  other  image  types.  Interpretation  of  image 
contents  apparently  relies  on  the  use  of  the  stereo 
Impression  so  that  a  pure  single-image  approach  is  not 
accepted  even  with  existing  DHMs. 


2.7  Past  Accuracy  Studies 

Numerous  authors  have  compared  single  radar  images 
to  maps,  thereby  employing  one  of  the  described 
computing  methods.  An  example  of  an  early  effort  was 
by  Leberl  (1972a),  as  shown  in  Figure  2.11.  Various 
interpolative  methods  were  used  to  describe  the  image 
deformation. 


FIGURE  2.11:  Example  of  a  data  set  for  an  early 
single  image  accuracy  elevation. 
British  real  aperture  radar  flight 
(1965)  over  Noord-Oost-Polder,  the 
Netherlands,  and  map  at  scale 
1:50  000  for  comparison. 


Table  I  is  a  review  of  some  existing  work.  In 
view  of  the  many  possible  mapping  methods,  of 
variations  in  the  image  types,  of  input  data  available 
on  the  terrain  etc.,  it  is  ovious  that  a  comparison  of 
results  is  nearly  impossible.  However,  one  conclusion 
can  be  drawn,  the  range  of  achievable  results  stretches 
from  the  geometric  ground  resolution  (Grade  et  at. 
1970)  to  some  large  figures.  Clearly,  though  these 
larger  figures  may  be  caused  by  unknown  topographic 
relief,  lack  of  ground  control  points  or  poor  internal 
geometry. 

Expectations  were  expressed  in  the  preparation 
phase  of  satellite  radar  that  it  could  result  in  higher 
image  geometry  performance  than  aircraft  radar  (Leberl, 
1976b).  SEASAT-SAR  was  the  first  data  to  potentially 
prove  that  expectation.  However,  the  overall  system 
engineering  for  SEASAT-SAR  did  not  optimize  geometric 
performance.  Table  II  presents  results  obtained  with 
SEASAT-SAR.  Hardcopy  images  on  film  either  from 
optical  or  from  digital  correlation  provided  the  same 
accuracy.  Time  and  range  referencing  was  not  of 
sufficient  quality  to  fully  capitalize  on  a  relatively 
stable  orbit  geometry  and  the  resulting  good  geometric 
potential  (Leberl  et  al.  1981).  Hithout  any  ground 
control  no  meaningful  positioning  of  ground  points  was 
feasible.  With  one  or  two  such  points  the  positioning 
error  is  ♦  300  to  400  m. 


TABLE  II 

SEASAT-SAR  geometric  accuracy.  Coordinate  errors  in  km  after  radargrammetric 
measurements  and  computation,  using  the  arctic  passes  over  Banks  and  Victoria  Islands 
taken  by  SEASAT  covering  3S  000  km2,  and  pass  3S1  over  Los  Angeles 


No.  of  control 
points 

Rigorous  radar¬ 
grammetry 

*  y 

Simple  scale  fit 

x  y 

Sea  arctic  images  of  Banks  Island,  35  000  sq  km 

1 

0.45 

0.49 

— 

— 

2  (at  both  ends) 

0.33 

026 

0.37 

0.37 

17 

0.17 

016 

- * 

— 

Image  over  Los  Angeles,  7  pts  per 

100  sq  km 

Optical  con. 

0.024 

0.017 

— 

— 

Digital  con. 

0.022 

0012 

— 

The  rigorous  use  of  digital  imagery  enables  one  to 
more  fully  exploit  the  data  geometrically  (Brown  et  al. 
1981).  Each  pixel  can  be  associated  with  a  precise 
sensor  events  time  and  a  slant  range  between  antenna 
and  object.  As  a  result  the  location  of  an  object  on 
the  Earth's  surface  can  be  computed  without  any  ground 
control  point  to  within  several  hundreds  of  meters. 
For  the  case  of  SEASAT ,  there  is  thus  little  difference 
between  the  use  of  rigorous  radargrammetry  and  coarse 
approximations,  except  for  the  total  lack  of  ground 
control  or  the  presence  of  merely  one  or  two  points. 
For  more  appropriate  mapping  applicability  a  satellite 
radar  would  have  to  be  equipped  with  precise  'time  and 
range  references  and  would  need  to  be  accompanied  by 
accurate  satellite  orbit  data  with  errors  smaller  than 
the  geometric  resolution  of  the  sensor. 


3.0  DESCRIPTION  OF  A  SYSTEM  FOR  DIGITAL  RADAR  IMAGE 
RECTIFICATION 


3 . 1  General 


Input  to  a  system  for  radar  image  rectification 
consists  of: 

-  a  digital  radar  image, 

-  sensor  interior  orientation  parameters, 

•  sensor  position  and  attitude  values, 

-  ground  control  points, 

-  a  raster  digital  elevation  model  (DEM). 

In  the  current  context  one  has  to  use  the  input  to 
first  create  a  mathematical  relationship  between  image 
and  map  and  then  to  actually  create  a  rectified  image. 
Many  alternatives  could  be  defined.  The  system 
developed  under  this  study  is  based  on  the  assumption 
that 

-  images  exist  of  mountainous  areas, 

-  identification  of  ground  control  point’s  is 

difficult. 

Therefore  the  system  is  based  on  radar  image  simulation 
to  support  the  generation  of  a  precise  correspondence 
between  radar  image  and  object  coordinates. 

Radar  image  simulation  with  DEMs  is  only 
meaningful  in  the  case  of  mountainous  terrain;  in  flat 
areas  a  simulation  as  used  here  would  be  redundant. 

A  simulation  system  SIMRISA  (Simulation  of  JRadar 
Images  using  an  ^mage  Space  Algorithm)  has  been 
developed  and  used  for  radar  stereo  analyses  (Domik, 
Leberl  and  Raggam,  1983).  It  is  extended  to  serve  for 
rectification  in  the  current  effort. 


3 . 2  Radar  Recti f ication  Based  On  Simulation 

The  procedure  used  for  the  rectification  of  radar 
images  of  mountainous  terrain  consist  of  the  following 
steps : 

(a)  Generation  of  a  digital  elevation  model,  either 
from  a  map  or  from  other  external  sources,  or 
from  radar  stereo  mapping; 

(b)  Verification  and  improvement,  or  computation  of, 
sensor  position  and  attitude  parameters  based  on 
identifiable  ground  control  points; 

(c)  Using  the  outputs  of  steps  (a)  and  (b),  i.e.  a 
DEM,  radar  imaging  parameters  and  a  flight  path, 
create  a  simulated  radar  image  to  resemble  the 
given  input  image; 

(d)  Correlate  the  real  and  the  simulated  radar  image, 
either  by  pointing  to  homologue  features  or  by 
automated  correlation,  thereby  creating  a  grid  of 
anchor  points  in  the  real  radar  image  to  relate 
to  the  simulated  one  and  consequently  to  the 
object  coordinate  system; 

(e)  Using  a  warping  function,  correct  the  real  image 
to  fit  over  the  simulated  image; 

(f)  Generate  a  geometrically  rectified  real  image  by 
assigning  the  real  image  gray  value  associated 
with  the  simulated  pixel  to  each  DEM  grid  cell; 

(g)  By  subtraction  of  the  real  image  gray  values 
from  the  simulated  image  one  obtains  a  radio- 
metrically  corrected  output  image. 

The  above  steps  imply  that  an  extensive  software  and 
radar  image  processing  system  is  available.  In  fact 
the  system  needs  to  consist  of  various  separate 
elements.  In  the  current  context  one  has: 

-  DEM  generation  by  the  Graz  Terrain  Model  (GTM) 
system  which  takes  given  map  contour  lines  to 
create  a  raster  DEM; 

-  System  for  stereo  mapping  with  radar  techni¬ 
ques  (SMART)  to  extract  contour  lines  or  raster 
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DEM  data  from  stereo  radar  imagery  on  a  photo- 
grammetric  analytical  plotting  machine  (Kern 
DSR-1  ) ; 

-  Radar  space  resection  programs  to  compute  the 
flight  parameters  from  a  given  input  image  and 
ground  control  points; 

-  Image  processing  programs  for  image  correlation 
display  and  coordinate  extraction; 

-  General  purpose  rectification  based  on  a  grid 
of  anchor  points  (RECTIF); 

-  Radar  image  simulation  (SIMRISA); 

-  Assignment  of  gray  values  to  DEM  grid  cells  based 
on  an  address  file  as  a  result  of  correlating 

the  real  with  the  simulated  image  (GEOREC). 

The  large  set  of  .'equired  software  for  this 
purpose  was  envolved  over  the  years  in  a  sequence  of 
individual  efforts.  The  following  will  describe  some 
of  the  more  relevant  components,  in  particular  the 
method  of  DEM  generation  from  stereo  radar  images  and 
the  radar  image  simulation. 


4  .0 


GENERATING  A  DEM  KITH  RADAR  STEREO  MAPPING 


4 . 1  Background 

Modern  photogramrae tri c  analytical  plotters  have 
opened  up  the  possibility  to  employ  other  types  of 
images  than  metric  aerial  photography  for 
stereo-plotting.  Instead  of  a  fixed  mechanical  or 
optical  analogon  for  the  projection  ray  one  has  a 
flexible  computer  to  relate  image  xy-coordinates  to  the 
cartesian  XYZ  object  or  model  system. 

Norvelle  (1972)  probably  was  the  first  to  actually 
program  an  analytical  plotter  for  stereo  radar.  The 
flexibility  of  earlier  analytical  plotters  was  hampered 
by  the  available  computers  and  the  use  of  low-level 
programming  languages.  Prior  to  Norvelle's  work 
stereo-radargrammetry  was  developed  in  the  context  of 
basic  studies,  e.g.  by  Rosenfield  (1968)  and  Gracie  et 
al.(1970)y  and  had  led  to  a  rather  unique  and  specific 
radar  stereo  mapping  instrument  (Yoritomo,  1972; 
Graham,  1972). 

The  current  generation  of  analytical  plotters 
employs  general  purpose  computers,  distributed 
processing  and  high  level  programming  languages. 
Therefore  the  effort  required  to  program  them  is  much 
reduced  when  compared  to  the  early  systems.  This  fact 
is  combined  with  an  increased  interest  in  radar  remote 
sensing  as  documented  by  satellite  and  other  radar 
imaging  experiments.  Therefore  it  appears  increasingly 
meaningful  to  study  radar  stereo  mapping. 

In  an  effort  to  demonstrate  the  flexibility  of 
current  commercial  analytical  photogrammetric 
instruments  and  in  order  to  explore  more  fully  the 
capabilities  and  limitations  of  stereo  radargrammetry  a 
procedure  and  software  system  for  the  analytical 
plotter  Kern  DSR-1  was  developed  to  set-up,  and  plot 
from,  radar  stereo  models.  This  is  a  first  and 
preliminary  report  on  this  development,  emphasizing  a 
broad  description  of  the  procedures  of  SMART  --  Stereo 
Mapping  with  Radar  Techniques . 1 ) 

The  chapter  will  discuss  the  procedure  and 
demonstrate  it  with  a  unique  data  set.  Common  radar 
image  acquisition  is  with  parallel  flight  lines.  It  is 
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an  added  complexity  to  use  intersection  flight  lines. 
As  will  be  shown  the  DSR-1  is  able  to  handle  this  type 
of  data  and  produce  a  digital  terrain  elevation  model 
of  an  accuracy  commensurate  with  the  radar  input. 


1) .  2  Overview  Of  System  Modules 

The  programs  for  stereo  radar  are  organized  in 
seperate  modules.  The  entire  process  of  model  set-up 
follows  radargrammetric  formulations  presented  earlier 
and  reviewed  by  Leberl  (1983).  Essentially  the 
procedure  requires  sufficient  ground  control  for  the 
transformation  of  radar  coordinates  to  the  object 
space.  The  individual  images  are  processed  first  to 
obtain  approximate  values  for  the  rigorous 
radargrammetric  solution.  In  the  terminology  of 
classical  photogrammetry  one  first  solves  resections  in 
space  for  the  two  images  to  then  continue  with  a  bundle 
adjustment  using  both  images  simultaneously. 

Should  no  ground  control  be  available  then  one  has 
to  work  with  given  or  assumed  parameters  of  the  flight 
path  and  sensor  attitude.  The  preliminary  resections 
would  be  skipped. 


(a)  Information  and  Initialization  Module?  This 
task  enables  the  input  and  manipulation  of  initial 
project  and  image  parameters.  It  serves  for  general 
information  for  the  operator,  such  as  for  example 
project  name,  available  images  for  the  project  or  image 
and  orbit  data  Information. 


( b )  Control  Point  Management  Module :  This  module 
is  for  the  input  and  manipulation  of  ground  control 
point  data.  The  ground  coordinates  are  in  an 
orthogonal  cartesian  coordinate  system,  usually 


This  development  was  done  outside  the  supported  part  of 
the  current  project  DA JA**5-83-C-00 1 8 .  It  is  reviewed  here 
due  to  its  interest  to  the  radar  rectification  and  for 
ease  of  reference;  see  also  Raggam  and  Leberl  (198*1). 
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referenced  to  an  origin  within  the  area  to  be  mapped 
and  the  Z-axls  pointing  along  the  local  vertical  of  the 
origin . 


(c) 

Orbit 

or 

Flight  Data 

Management 

Module : 

This 

is 

for 

input 

an 

d  manipula 

tion  of  sensor  position 

mea 

sureme 

nts . 

The 

system  ass 

umes  the  orbit  or 

flight 

to 

be  re 

presen 

ted 

with  time 

polynomials 

so  that 

sensor 

pos 

itions 

s  may 

be 

expressed 

as  follows: 

where  t  is  the  time.  Approximations  for  the  orbit 
coefficient  vectors  (a,  t>,  c,  ..)  will  be  determined  in 
the  orbit  data  management  module.  Input  are  the 
position  measurements,  output  is  an  array  of  polynomial 
coefficients.  Joint  low-order  polynomials  will  be  used 
instead  of  higher  order  polynomials  if  the  need  exists. 


( d )  Single  Radar  Image  Processing:  This  module 
represents  the  equivalent  of  a  resection  in  space  to 
compute  the  orientation  parameters  for  a  single  radar 
image  with  inner  and  exterior  orientation.  Plotting 
with  a  single  radar  image  may  be  done  on  either  a 
spherical  surface  of  chosen  radius  or  in  a 
three-dimensional  X,  Y,  Z  coordinate  system  with 
known  constant  Z. 


( e )  Stereo  Radar  Image  Processing:  This  module 
Includes  the  actual  set-up  of  a  radar  stereo  model.  It 
follows  the  use  of  the  single  radar  image  module  to 
determine  approximations  of  the  exterior  orientation 
parameters  of  the  two  overlapping  images.  After 
elimination  of  y-parallaxes  this  module  serves  to  plot 
contour-lines  or  planimetric  features.  Results  are 
either  directly  plotted  on  an  xy  pen  plotter  or  are 
entered  into  a  digital  data  base  for  further 
processing.  One  example  that  is  currently  operational 
is  for  digital  terrain  elevation  models  with  the  Graz 
Terrain  Model  (GTM)  program  system  (Oswald,  Raetzsch, 
in  print). 
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^•3  Radar  Stereo  Model  Set-up 


As  previously  described  the  radar  stereo  model 
set-up  in  SMART  is  realized  in  two  steps  consisting  of 
the  computation  of  the  elements  of  inner  orientation  of 
the  radar  image  pair,  followed  by  an  exterior 
orientation  with  a  radar  bundle  adjustment. 


Inner  Orientation :  The  inner  orientation  requires 
the  establishment  of  a  relationship  between  plate  or 
image  x,  y  coordinates  and  the  physical  radar 
measurements  of  time  t  and  slant  range  r  (Leberl, 
1983)*  The  system  uses  a  so-called  'range  reference 
line'  at  the  near  range  edge  of  the  image;  it  would 
correspond  to  the  start  of  the  sweep  on  the  image 
recorder.  It  should  be  defined  by  some  distinct  tick 
marks  (fiducials).  In  the  event  that  no  such  marks 
exist  the  operator  has  to  create  artificial  marks  and 
determine  the  inner  orientation  in  a  process  of 
self-calibration  together  with  the  exterior 
orientation . 


Exterior  Orientation :  Measurement  of  homologue 

orientation  points  can  first  of  all  be  as  in  a 
comparator  by  removing  manually  the  x-  and 
y-parallaxes.  Secondly  after  a  preliminary  resection 
with  each  image  one  is  in  a  nearly  parallax-free  stereo 
model  and  observation  of  homologue  points  is  made  more 
convenient.  In  contrast  to  analog  stereo 
photogrammetry  the  analytical  plotter  does  not  need  a 
relative  orientation  task  seperate  from  absolute 
orientation.  The  parameters  of  relative  and  absolute 
orientation  are  found  simultaneously  with  a  so-called 
bundle  adjustment.  In  our  system  this  approach  is  the 
one  used  for  radar  images.  The  basic  equations  for  the 
radar  stereo  adjustment  consist  of  two  types: 

(a)  The  squint  angle  condition: 

5-  (x?-  z.)  -  •  I  i  Mc-xf  *4 
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(b)  The  range  condition: 
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for  slant  range  presentations 
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for  ground  range  presentations. 

Each  measured  point  gives  rise  to  two  equations  of 
type  (  ^  .  1  )  and  ( *» .  2  )  .  In  addition  each  pair  of 
homologue  image  points  produces  one  additional 
condition.  This  is  called  the  cocircularity  condition, 
which  defines  the  location  of  an  imaged  point  in  object 
space  by  the  intersection  of  two  circles,  each  defined 
by  two  equations  of  type  { H . 1 )  and  (M.2).  The 
equations  are  non-linear.  Linearized  forms  are  used  in 
an  iterative  solution. 


Observations  for  the  computations  are  sensor 
position  measurements,  image  coordinate  measurements  of 
ground  control  points  in  one  or  in  both  images  and 
pairs  of  homologue  image  coordinate  measurements  of 
stereo  model  orientation  points.  Unknowns  to  be 
determined  within  the  adjustment  are: 

-  the  coefficients  of  the  orbit  time  polynomials; 

-  the  radar  imaging  parameters  (inner  orienta¬ 
tion)  to  convert  radar  image  coordinates  to 
range  and  time; 

-  the  parameters  of  a  correction  polynomial  to 
describe  radar  image  deformations; 

-  a  value  for  the  squint  angle. 


The  solution  for  the  unknowns  in  each  iteration  of 
the  non-linear  equation  system  is  obtained  in  the 
instrument's  host  processor  in  a  least  squares 
adjustment  using  the  method  of  conjugated  gradients 
(Schwarz,  1970).  This  solution  method  is  iterative  and 
may  have  some  advantages  for  larger  equation  systems 
when  a  small  computer  must  be  used  and  both  computing 
times  and  storage  requirements  are  limited. 
Furthermore  there  is  evidence  that  the  method  is  more 
robust  for  this  special  problem  than  a  strict  Gauss 
solution  of  normal  equations. 

Depending  on  the  quality  of  the  approximate  values 
the  adjustment  process  must  be  repeated  iteratively 
until  results  satisfy  specified  termination  criteria. 
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-  th<  location  of  the  pixel  in  the  ir.ape. 

The  gray  value  is  a  function  of  incidence  angle, 
wavelength,  polarisation  and  reflecting  capability  of 
the  corresponding  targets.  With  fixed  wavelength  and 
polarisa  tion  the  brightness  varies  only  with  incidence 
angle  and  reflection  properties. 

The  location  of  the  pixel  in  the  image  is  defined 
by  the  time  delay  of  the  radar  pulse  to  and  from  the 
target.  For  a  perfectly  flat  ground  the  time  delay 
increases  proportional  to  the  distance  of  nadir  and 
target.  Irregularities  on  the  surface  to  be  imaged 
affect  the  time  delay  and  create  distortions  in  the 
image . 


An  extreme  effect  of  relief  displacement  occurs 
with  lay-over;  the  top  of  a  traget  to  be  imaged  is 
closer  to  the  radar  sensor  than  the  bottom  and  is 
therefore  recorded  sooner.  However,  lay-over  is  also 
defined  as  the  simultaneous  mapping  of  two  different 
ground  targets  into  one  resolution  cell.  This  is  the 
case  when  the  slant  range  of  two  objects  is  the  same. 
Lay-over  effects  the  brightness  as  the  return  power 
from  two  or  more  targets  is  added  up.  The  simulation 
process  can  correctly  consider  the  lay-over  case  by 
attaching  more  than  one  position  address  to  an  image 
gray  value. 

After  the  simulation  and  registration  each  digital 
DN  number  of  the  corrected  radar  image  is  connected  to 
one  or  more  DEM  addresses.  The  procedure  GEOREC 
(GEOmetric  RECtif ication )  performs  the  final  step  of 
the  radar  image  rectification  making  use  of  the  data 
structures  described  in  chapter  5.3. 

Input  to  this  last  step  are 

-  grey  values  of  the  real  image  registered 
to  the  simulation 

-  address  file  (containing  DEM  addresses) 

-  relation  between  grey  values  and  address 
file 
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6.0  OTHER  AUXILIARY  SYSTEM  ELEMENTS 


6 . 1  Re  r,  1st  ration 

The  image  registration  is  based  on  a  program 
RECTIF  by  Diarra  (1982);  it  is  a  modular  software 
system  for  geometric  correction  of  an  image.  Input  is 
a  set  of  honologue  pairs  of  points  whose  coordinates 
are  known  in  the  real  image  and  in  the  simulation. 
These  are  used  to  compute  correction  polynomials  to 
describe  the  trend  of  the  image  deformations. 
Residuals  at  control  points  serve  then  to  compute 
refined  corrections  with  some  higher  order 
interpolative  method.  Work  is  thus  via  "warping"  or 
"rubber-sheeting".  Polynomials  are  limited  to  degree  3 
or  less,  so  that  they  are  both  stable  and 
overdetermined.  Residuals  are  extracted  at  each  anchor 
point,  and  serve  to  modify  the  original  polynomials. 

These  correction  polynomials  are  used  for 
resampling  based  on  a  grid  of  anchor  points;  each  set 
of  object  coordinates  obtains  image-positions.  Various 
methods  may  then  be  applied  to  determine  gray  values  at 
those  positions:  interpolation  followed  by  nearest 
neighbor  assignment,  bilinear  interpolation,  and 
bicubic  interpolation,  are  the  current  options. 
Control  points  will  usually  be  defined  by  hand,  or  by 
digital  correction. 

The  registration  results  in  a  rectified  image;  in 
the  case  of  data  actually  studied  remaining  distortions 
are  of  an  average  of  3  pixels  or  less,  assuming  an 
appropriate  distribution  of  control  points. 

This  results  in  a  real  radar  image  superimposed 
onto  the  simulated  data. 


6.2  Removing  Geometric  Radar  Distortions 

Any  gray  value  occuring  in  the  digital  radar  image 
is  characterized  by: 

-  the  gray  value  number  (brightness)  of  the  pixel; 
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5 . 4  Software  Implementation 

The  product  of  the  software  implementation  to 
perform  the  simulation  is  named  SIMRISA  (SIMulation  of 
Radar  images  using  an  Image  Space  Algorithm). 

Emphasis  was  put  on  modularity  of  the  system  (see 
description  of  system  modules,  Fig.  5.3)»  The  program 
can  be  easily  changed,  e.g.  additional  thematic  input 
data  sets  could  allow  for  use  of  different  baekscatter 
curves  or  use  of  different  look-up  tables  within  one 
DEM  area. 


In  areas  other  than  shadows  the  inape  g,ray  value 
is  a  function  of  the  incidence  angle  for  the  incoming 
radiation.  This  is  an  angle  between  the  line  of  sight 
antenna-object  and  the  local  normal  onto  the  surface. 

Once  this  angle  is  computed  for  a  pixel  one  has 
the  choice  of  selecting  a  back-scatter  curve  according 
to  thejtheroatlc  charactaristics  of  the  material  and  of 
the  ylngle  4*  .  Standard  backscatter  function  are  those 
of  Hagfors  (196*1),  Muhleman  (196*1)  or  a 
cosine- function. 

Simultaneous  illumination  of  various  terrain  parts 
leads  to  image  lay-overs:  the  projection  circle 
intersects  more  than  one  point  on  the  ground.  The  gray 
value  is  the  sum  of  the  individual  reflection  values. 


5. 3  Data  Structure 

The  internal  data  structures  are  constructed  to 
make  available  a  fast  access  to  the  address  file 
elements  in  relation  to  the  image  elements.  Each 
azimuth  line  in  the  image  represents  a  logical  block 
and  corresponds  to  a  unique  logical  block  of  the 
address  file.  Memory  allocation  for  each  logical  block 
will  be  a  minimum  of  n*e  bytes,  where  n  is  the  number 
of  range  columns  and  e  the  number  of  bytes  used  for 
storage  of  one  gray  value.  Address  file  B  allocates 
(at  minimum  again)  h*f  bytes,  with  h  being  the  number 
of  addresses  belonging  to  one  azimuth  line  and  f  the 
number  of  bytes  used  for  one  address  definition.  The 
relation  is  realized  by  an  additional  access-matrix 
(pointer  file),  which  allows  a  fast  association  between 
elements  of  the  two  different  files.  The  additional 
storage  for  the  relation  matrix  will  be  h*2  bytes  per 
logical  block. 
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The  sensor  position  and  velocity  vector,  derived 
from  the  i macing  time,  define  a  cartesian  coordinate 
system:  The  sensor  position  is  the  origin  of  this 
"sensor  coordinate  system",  the  axes  are  defined  by  the 
velocity  vector  and  the  sensor  position  as  expressed  in 
a  geocentric  coordinate  system  (Eq.  2.8).  The  set  up 
of  this  system  allows  a  clear  presentation  of  the 
geometry  created  by  any  input  coordinate  (i,j). 

The  origin  of  the  sensor  coordinate  system  is  the 
center  of  a  sphere  with  radius  r  (range  sphere).  In 
general,  a  squint  angle  'V  t  4  is  assumed  for  imaging. 
Thus  a  cone  is  defined  by  the  sphere  and  its  axis  along 
the  velocity  vector:  if  the  squint  angle  TT  is  equal  to 
4  Imaging  is  in  a  plane  perpendicular  to  the  velocity 
vector.  The  intersection  of  sphere  and  cone,  or, 
sphere  and  plane,  respectively,  defines  a  circle 
(intersection  circle).  The  remaining  task  is  now  to 
intersect  the  circle  with  the  DEM. 

At  least  one  intersection  point  but  possibly 
several  ones  will  be  found.  The  occurrence  of  more 
than  one  solution  is  called  "lay-over":  different 
terrain  targets  were  illuminated  simultaineously .  The 
address  of  each  Intersection  point  as  well  as 
corresponding  look  angle,  surface  inclinations  and 
additional  informations  are  stored  for  further 
processing. 


5.2  Radiometric  Model  Of  Imaging 

The  definition  of  an  image  gray  value  in  a  given 
pixel  location  is  denoted  as  "radiometric  model".  One 
has  to  define  areas  of  shadow  (black),  specular 
reflection,  multiple  reflection  in  lay-overs  etc. 

Black  image  areas  are  mostly  shadows  or  water 
bodies.  In  order  to  find  those  pixels  that  indeed  are 
shadow  areas  one  has  to  sort  the  look  angles  going  from 
near  range  to  far  range  along  an  imaged  object  line. 
Shadow  areas  begin  where  an  increasing  range  r  is 
associated  with  a  reduced  look-angle  off  nadir. 
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5.0  RADAR  IMAGE  SIMULATION 

Approaches  to  image  simulation  can  be  categorised 
according  to  the  domain  of  the  transformation  function. 
Object  space  algorithms  map  each  object  point  (DEM 
point)  into  a  two  dimensional  field  (image)  whereas 
image  space  algorithms  start  from  equidistant  image 
coordinates  (i,j)  of  the  output  image  and  assign  gray 
value  accordingly  to  it.  Object  space  algorithms  are 
usually  faster  in  their  application  than  image  space 
algorithms.  But  there  is  a  need  for  interpolation  to 
convert  the  unequally  spaced  image  addresses  to 
equidistant  ones. 

The  development  of  the  radar  simulation  was  part 
of  another  project  (see  Domik,  Leberl  and  Raggam, 
1983).  It  is  of  great  relevance  in  the  current  context 
and  is  discussed  in  the  following  for  case  of 
reference . 

An  image  space  algorithm  was  chosen  to  be 
implememted  because  it  provides  DEM  addresses  at  exact 
image  locations.  Fig.  5.1  and  5.2  show  the  basic  idea 
of  object  vs.  image  space  algorithms. 

Gray  values  are  attached  to  the  imagr  coordinates 
dependent  on  the  corresponding  incidence  angle  on  the 
ground. 


5 . 1  The  Geometric  Imaging  Model 

The  image  coordinates  (i,j)  serve  as  start  values 
for  the  simulation  program.  They  are  used  to  extract 
imaging  time  t  and  slant  range  r  or  ground  range  g  by 
Eq .  '  s  5.1* 


t  = 

+  *0 

EQ. ( 5 . la ) 
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+  sd 

EQ. (5.1b) 
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l< .  5  Software  Implementation 

The  software  package  is  currently  installed  at  a 
Kern  DSR-1  analytical  plotter  under  the  name  SMART 
(Stereo  Mapping  with  Radar  Techniques). 

The  DSR-1  standard  configuration  consists  of  three 
independent  digital  processing  units,  which  are  denoted 
as  PI,  P2  and  P3  (Chapuis,  1980).  Processor  PI  is  the 
host  computer.  It  serves  for  development  and  execution 
of  application  programs  and  performs  data  transfer 
between  the  other  processors.  The  computation  of 
set-up  parameters  for  the  radar  stereo  model  is  done  in 
PI.  Processor  P2  receives  the  radar  image  orientation 
parameters  from  PI  and  uses  them  to  convert  in  real 
time  model  coordinates  X,Y,Z  to  radar  image  coordinates 
r,,t,,r,,,t"  and  onto  DSR-1  plate  coordinates  x*,y'  and 
x",y".  P2  also  transfers  image,  model  or  object 
coordinates  to  processor  PI  if  requested.  Processor  P3 
is  an  added  convenience  for  the  operator  and  controls 
an  operator  control  panel  which  may  be  used  for 
communication  between  operator  and  DSR-1. 

The  Kern  DSR-1  therefore  uses  three  computers  as 
opposed  to  other  plotters  which  may  employ  only  one. 
What  seems  to  be  an  added  complexity  in  fact  is  an 
advantage:  the  real-time  operations  are  performed  by  a 
seperate  processor  P2  with  a  clearly  defined  software 
interface  for  a  user  working  with  the  host  PI.  The 
user  will  therefore  work  in  a  familiar,  general  purpose 
multi-user  environment  with  an  unmodified  operating 
system.  The  critical  real  time  operations  are  singled 
out . 


mmkxkkxkkxxxkxxkkxkkkxkkxxkkkkkxkkxx 

*  CEPHALONIA  AND  ITHAKA  x 

*  data  generated  at  kern  dsr-i  k 
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FIGURE  *».3:  Digital  elevation  model  of  the  greek 
islands  Cephalonia  and  Ithaka  in 
axonometric  view.  Data  digitized  at 
KERN  DSR-1  stereoplotter  using 
program  SMART. 
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^  *  **  Example  Of  Radar  Stereo  Plotting  With 

Satellite  Da ta 

A  Space  Shuttle  Imaging  Radar  SIR-A  stereo  image 
pair  of  the  greek  islands  Cephalonia  and  Ithaka  was 
used  to  illustrate  the  procedure  of  extracting  a 
digital  height  model  (Fig.  it .  1  > .  A  more  detailed 
study  of  the  stereoscopic  computations  from  such  data 
is  presented  by  Kobrick  et  al .( submitted ) .  Viewing 
limitations  were  discussed  by  Domik  et  al.(1983). 
After  model  set-up  the  topographic  elevations  were 
digitized  in  the  object  coordinate  system  using 
profiles  and  ridge  lines;  the  measurements  were 
entered  into  the  terrain  model  system  GTM.  The  GTM 
program  enables  one  to  generate  contour  line  plots  at 
chosen  interval,  to  obtain  axonometric  views  of  the 

terrain  data  and  various  other  products. 

Figure  ^.2  shows  an  axonometric  view  of  the 

digital  elevation  model  (DEM)  of  the  islands  extracted 
from  a  map  1  :  200  000,  Figure  M . 3  presents  the  DEM 
created  on  the  Kern  DSR-1  analytical  stereoplotter.  A 
comparison  of  the  two  DEM's  is  possible  with  the 
GTM-system  and  reveals  that  the  root  mean  square 
differences  over  all  DEM  points  of  the  islands  amount 
to  +  98  meters. 


According  to  the  data  structure  description  the 
logical  blocks  of  grey  value  and  address  file 
correspond  sequentially  to  each  other.  This  allows  to 
process  each  azimuth  line  (corresponding  to  one  logical 
block)  independently.  Each  gray  value  is  now  placed  at 
the  DEM  addresses  calculated  during  the  simulation. 
Thus  the  geometric  radar  distortions  are  eliminated 
again  and  the  resulting  ortho  image  displays  distances 
instead  of  (slant  or  ground)  ranges. 

Usually  the  DEM-addresses  are  not  exact  grid 
positions,  so  that  some  kind  of  interpolation  has  to 
take  place. 


7.0  EXAMPLES  OF  RADAR  IMAGE  RECTIFICATION' 


7 . 1  Oetztal-Data 

A  radar  image  was  obtained  during  the  European 
SAR-580  campaign  of  the  Oetztal  in  Tyrol  (Rott,  1983). 
The  sensor  for  this  mission  was  carried  aboard  a 
Convair-580  aircraft  with  a  dual  band,  dual 
polarisation  SAR  which  could  be  operated  in  X  and  C 
band  or  X  and  L  band.  Data  were  recorded  by  using 
either  an  optical  or  a  digital  recorder: 

The  aim  of  the  Tyrol-part  of  the  SAR-580 
experiment  were  investlgatlors  in  the  applications  of 
radar  for  mapping  snow  and  glaciers  in  mountainous 
regions.  The  sensor  was  operated  in  X  and  C  band  (HH 
polarisation),  covering  a  ground  swath  of  13  km  length. 
Fig.  7.1  is  a  presentation  of  the  SAR  image,  X  band, 
recorded  digitally  and  processed  at  the  German 
Aerospace  Research  Establishment  (DFVLR).  The 
presentation  is  in  ground  range;  the  original 
resolution  of  3m  was  changed  to  a  coarser  resolution 
of  24  m  for  presentation  in  512  x  512  pixels. 

A  DEM  was  derived  from  a  map  1  :  50  000  by 

scanning  the  contour  and  drain/ridge  lines.  The  imaged 
part  is  only  a  subarea  of  the  DEM,  which  extends  from 
1 0°98 '  to  10°59 *  East  longitude  and  from  H6°95* 
to  A6°52*  North  latitude.  Fig.  7. 2a  shows  an 

axonometric  presentation  of  the  height  model.  The 
mountain  spreading  from  the  south-western  towards  the 
north-eastern  corner  can  be  recognized  in  the  radar 
image  (upper  half). 
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7 . 2  Processing  Of  The  Oetztal  Data 

The  simulation  was  carried  out  by  the  simulator 
according  to  chapter  5.  Since  flight  navigation  was 
eroneous  during  the  SAR-580  flight,  one  has  to  rely  on 
control  points  to  approximately  derive  the  flight  and 
imaging  parameters  (Fig.  7.3).  The  maximal  height  in 
this  area  is  3615  m,  whereas  the  minimal  height  is  85 
m.  Hagfors  law  was  used  to  model  the  backscatter 
depending  on  the  incidence  angle.  No  additional  input 
was  used  to  distinguish  between  different  thematic 
features  in  the  image. 

The  simulation  results  in  a  gray  value  image  with 
an  additional  address  file  and  relation  matrix.  The 
gray  values  assigned  by  the  simulator  are  merely 
functions  of  the  incidence  angle  and  do  not  consider 
any  properties  of  the  imaged  targets  (Fig.  7. 4). 

The  geometry  created  by  the  simulation  matches  the 
real  image  closely.  The  radar  distortions  due  to  radar 
properities  like  forshortening  and  shadowing  are 
correctly  implemented  by  the  simulation. 

However,  sensor  dependent  errors  will  be  present 
due  to  the  Inaccurate  flight  and  imaging  records.  Thus 
there  is  a  need  for  geometric  image  correction.  The 
registration  described  in  chapter  6  creates  a  modified 
radar  image  to  fit  with  the  simulated  one.  The 
simulated  image  serves  to  support  control  point 
definition. 

The  corrections  applied  to  the  real  radar  image 
are  shown  in  Fig.  7.5*  Corresponding  features  in  both 
real  and  simulated  image  were  used  as  control  points  in 
the  registration.  Their  residual  vectors  are  plotted 
into  the  deformation  grid  with  a  10  x  10  pixel  size. 
The  DN  values  of  the  new  SAR  image  (Fig.  7.6)  may  now 
replace  the  synthetic  gray  values  and  correspond  to  the 
additional  output  files  of  the  simulation.  Thus  the 
connection  between  radar  image  and  DEM  addresses  is 
made  and  the  rectification  can  be  carried  out  by  the 
appropriate  procedure. 

The  radar  ortho  image  is  shown  in  Fig.  7.7. 
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The  sinulation  displays  the  gray 
function  of  the  incidence  angle,  wh 
values  in  the  real  inage  vary  with 
incidence  angle  and  thematic  c 
Substracting  the  simulated  gray  values 
data  leaves  a  radiometric  corrected  ima 
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lighter  gray  to  white.  Dark  gray  tones 
strong  reflection  in  the  simulation  as 
real  image.  Especially  near  the  edges  o 
dark  tones  are  largely  insignificant. 
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7 . 3  Cepha Ionia  Data 


On  Nov.  12,  1981,  the  SIR-A  (Shuttle  Imaging 

Radar-A)  was  launched  aboard  NASA's  space  shuttle. 
SIR-A  is  a  left-side-looking  SAR  operated  at  L-band 
frequency  with  an  elevation  angle  of  approximately  50 
off  nadir  to  the  center  of  the  imaging  swath. 
SIR-A  data  were  optically  recorded. 

Two  images  obtained  in  Europe  from  Data  Take  32-33 
and  Data  Take  37A,  resp. ,  were  of  special  interest  for 
stereo  investigations  (see  Fig.  4.1).  Though  SIR-A 
was  not  specially  designed  for  stereo-investigations, 
several  data  takes  from  crossing  orbits  proved  to 
deliver  fairly  well  stereo.  At  the  crossing  of  the 
above  mentioned  Data  Takes  two  Greek  islands, 
Cephalonia  and  Ithaca,  were  imaged.  The  intersection 
angle  of  the  imaged  swaths  is  approximatley  34°. 

For  our  investigations  the  images  presented  on 
film  were  scanned  by  an  optronics  device  to  receive  the 
data  in  digital  form.  The  digital  data  was  processed 
to  receive  a  resolution  of  100m. 

The  analog  images  were  used  to  create  a 
stereo-derived  elevation  model  as  explained  in  chapter 
4.  At  the  same  time  an  elevation  model  was  genereated 
by  digitizing  the  contour-lines  at  a  200  m  interval  of 
a  1  s  200,000  map  in  UTM  projection.  Axonometric 
presentations  of  both  files  are  presented  in  Fig.  4.2 
and  Fig.  4.3. 


7.4  Processing  Of  Cephalonia  Data 


The  process  is  analog  to  the  Oetztal  data 
processing,  using  the  modules  of  the  radar  image 
processing  system  described  in  chapters  4  to  6. 


Flight  and  imaging  information  to  be  entered  into 
the  simulation  were  gained  from  (Cimino,  198?). 
Additional  control  data  to  correct  the  parameters  were 
not  needed  in  this  case.  Through  the  simulator  a 
synthetic  radar  image  was  created  (Fig.  7.9)  from  the 
stereo-derived  DEM  together  with  an  additional  address 
file.  The  registration  results  in  Fig.  7.10  by 
removing  the  differences  between  real  and  simulated 
image.  Fig.  7.11  shows  the  ortho  image  after 
elimination  of  the  range  distortion. 
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CONCLUSION  AND  OUTLOOK 


We  have  described  an  elaborate  software  system  for 
the  rectification  of  digital  radar  images. 
Rectification  is  based  on  digital  elevation  models  from 
maps  or  from  radar  stereo  models,  on  control  points  and 
on  given  data  of  the  sensor  and  flight  path. 

The  core  of  the  procedure  is  an  image  simulation 
to  support  the  identification  of  homologous  features  in 
both  the  image  and  the  DEM.  Rectification  can  be 
geometric  and/or  radiometric. 

The  entire  system  is  experimental:  the  size  of 
manageable  data  sets  is  limited,  required  computing 
resources  are  significant.  In  order  to  improve  the 
current  system  one  has  to  speed  up  the  processes  and  to 
gain  more  experiences  with  both  the  simulation  and  the 
rectification.  Simulations  must  improve  to  consider 
thematic  variations,  and  to  more  distinctly  be 
sensitive  to  geometric  resolution. 

A  question  of  great  concern  is  the  matching  of 
images  from  different  sources,  such  as  Landsat-MSS,  -TM 
and  radar.  It  is  thus  clear  that  great  efforts  will 
still  be  required  to  expand  current  capabilities  so 
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